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We currently perform calculations with an improved staggered fermion action (p4fat3). We use a 
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in the low temperature regime. We also comment on the quark mass dependence of the equation 
of state at high temperature. 
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1. Introduction 

One major topic in thermodynamics of QCD and, especially, lattice QCD is the calculation the 
Equation of State (EoS). While there has been great progress in the recent years still several issues 
remain unresolved. One of them is the observed deviation in the EoS of the lattice computation 
compared to what one expects from a Hadron Resonance Gas model (HRG). As can be seen on 
the right hand side of Fig. [I| the HRG description shows deviations from the lattice QCD data [|l|], 
especially in the low temperature regime. On the contrary, because the HRG model is in quite good 
agreement with experiment in the low temperature region |^, one would expect that the HRG 
model also gives a good description for the low temperature region of the lattice QCD equation of 
state. A possible source for the discrepancy, aside from lattice discretization errors which still may 
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Figure 1: Trace anomaly from lattice QCD at various Nt for = O.lm, (left). Deviations are observed 
between HRG and lattice QCD for the low temperature regime (right). 

play a role on the = S lattices, may be due to the fact that previous lattice studies of the EoS had 
been performed with light quark mass values which are about a factor two larger than in nature. 
This corresponds to a light pseudo-scalar mass of about 220MeV. In order to quantify the effects of 
too large light quark masses we present new calculations for the lattice QCD equation of state at a 
smaller light quark mass niq = 0.05ms compared to our earlier calculation with = 0. Irris where 
nis is the strange quark mass. 

The new data for niq = O-OSm^ also is used to investigate the properties of the chiral phase 
transition. In particular, we will present results on the scaling with the light quark mass of the 
chiral condensate and the chiral susceptibility which is relevant for identifying the critical point of 
the chiral phase transition. 

The paper is organized as follows. In Sec. ^ we give an overview of numerical details and 
outline the calculation of the EoS on the lattice. In Sec. ^ the calculation of the potential and the 
change in the scale for the = 0.05mi case is discussed. The results for the EoS with rrig = 0.05ms 
are presented in Sec. ^ and results for the chiral condensate and susceptibility are shown in Sec. ^ 
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2. Numerical Details and Setup 

The computation of the QCD equation of state on the lattice requires large resources. Our 
calculations have been performed on IBM BlueGene/L and QCDOC supercomputers at Lawrence 
Livermore National Laboratory, the New York Center for Computational Sciences (NYCCS), and 
Jiilich Supercomputing Centre. For our calculations we use the p4fat3 action and the RHMC algo- 
rithm [|l|]. The computational details are summarized in Tab. |l| and Ref. |^. Further information 
can be found in Refs. [|, |, i]. 



Masses: 


■ 


= O.lni.5 


m^ = 0.00081 - 


0.00370 




Mg ■ 


= 0.05m, 


m^ = 0.00120 - 


0.00145 


Pion Masses: 


mq ■ 


= O.lm, 


nijt w 220 MeV 






nig : 


= 0.05m, 


m;t « 160 MeV 




Volume (r / 0): 


nig ■■ 


= 0.1m. 


32^ X 8 






nig : 


= 0.05m, 


32^ X 8 




Volume {T = 0): 


nig : 


= 0.1m, 


324 






nig : 


= 0.05m, 


324 




Stats. {T / 0): 


nig : 


= 0.1m. 


# 8,000-37,000 






nig : 


= 0.05m, 


# 5,000-20,000 




Stats. (T = 0): 


nig : 


= 0.1m. 


#2,000 - 6,000 






nig : 


= 0.05m, 


# 1,500-2,000 





Table 1: Overview of the simulation details. 



In the following we give a brief sketch of the calculation of the EoS on the lattice. The basic 
quantity is the trace anomaly 0'^'^ =e — 3p with the energy-momentum tensor Q''^, energy density 
e, and pressure p. The trace anomaly can be computed on the lattice in terms of the action density, 
and the strange and light quark chiral condensates. We are only interested in the thermal part of the 
trace anomaly and subtract the zero temperature part, 

&^^=e-3p={e- 3p)t -{e- 3p)t=o- (2.1) 

From this expression we can determine the pressure through the thermodynamic relation jj'' = 
Td{p/T^)/dT which is integrated to yield 

P{T) PiTo)_ [\fl_^,,^fy (2.2) 



r4 Jt, f5 

Here Tq is chosen to be in the deep hadronic region where p{To) is already exponentially small. 
The energy density e then is determined given and p. For further details see Ref. 



3. Potential 

To account for the quark mass dependence of the lattice cut-off, which corresponds to a shift 
in the temperature scale when we lower the light quark mass, we have calculated the potential for 
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Figure 2: Potential V{r) vs. improved distance rj at fi — 3.49 for = O.OSm., and niq — Q.lnis. Both 
quantities are given in dimensionless units using the Sommer scale parameter ro. 

several values of jS. For rriq = O.lniv the scale was calculated on 16^ x 32 lattices, see Ref. In 
the case of niq = 0.05ms we computed the potential on 32^ lattices. We have set ro = 0.469 fm. The 
results and some simulation details are summarized in Tab. ^ where we also compare lattice scales 
obtained for nig = 0.05m.s to niq = O.lm,. For illustration we plotted the potential for j3 = 3.49 in 
Fig. ^ both for mq = O.OSm^ and mq = O.lm^. 

From the preliminary results in Tab. ^ we conclude that within error bars the scale determined 
for rriq = 0.05ms is consistent with the older results for rriq = O-lm,. We are planning to increase 
statistics and include more j8 values in the future. For the moment we will make use of the global 
mass fit 

fl'm,=o.05m, = exp(-2A Am^) am^=o.lm,, , (3.1) 

in order to determine the scale. The value for A is found to be A = 1.40(3). In this fit more data 
points are involved than in our current analysis for rriq = O-OSnig. Using this fit we find a shift in 



p 


3.49 


3.51 


3.53 


3.54 




0.1455 


0.1370 




0.1272 


'^m,=0.05ms [fm] 


0.1435 


0.1385 


0.1306 




ro/o m, =0.1m.. 


3.223(23) 


3.423(30) 




3.689(37) 


ro/a m, = 0.05/11,, 


3.2675(90) 


3.386(20) 


3.592(11) 




rm,=0.1m, [MeV] 


169.5(1.2) 


180.0(1.6) 




194.0(1.9) 


Tmg=Q.05ms [MeV] 


171.85(47) 


178.1(1.1) 


188.90(60) 




^mij=0.\mi, 


330 


140 




300 


^m,=0.05mj 


200 


200 


170 





Table 2: Results and overview for the calculation of the potential. The last two rows give the number of 
configurations we have used in our analysis. 
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temperature of about 0.5 MeV. 



4. Mass Dependence of the QCD Equation of State 

In Fig. ^ we plot our result for the trace anomaly for ruq = 0.05mj. We compare to the niq = 
OAnis case from Ref. [[I], |3|] as well as to the HRG model result. The curves in the plot correspond 
to fits with a quadratic ansatz in T which fits the data quite well. Note that for our calculation of 
the trace anomaly for the HRG all resonances were chosen to be at their physical mass. We also 
note that for the = O.OSm^ case the Ughtest pseudo-scalar mass^ is about 160 MeV which is very 
close to its physical value. The shift of about 0.5 MeV in the scale, as discussed in Sec. ^, for the 
rriq = 0.05m.v case has been taken into account. 
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Figure 3: Trace anomaly for mq = 0.05ms compared to the HRG and nig — 0. Inis results. The curves are 
fits to the corresponding data made to highlight the cut-off and quark mass dependence of the various data 
sets. 

From Fig. |3| we observe in the transition region a shift towards lower temperature for {e — 
3p)/T^ for mq = O.OSm^ of about 4 MeV compared to niq = O.hris (with A^^ = 8)- For the larger 
temperature regime we find that for the smaller quark mass there is quite good agreement with the 
HRG model. However, for the lower temperature region we observe deviations from the HRG in 
that case. We will increase the statistics and add more data points in the lower temperature region 
in the future in order to clarify the situation. 

'Note that within the staggered formalism only one light pseudo-scalar exists at finite values of the lattice cut-off. 
Full flavor symmetry is only recovered in the continuum limit. 
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5. Chiral Condensate 



In this section we present new results on the chiral condensate and its susceptibility. We con- 
sider a combination of the light and strange quark condensate to eliminate quadratic divergencies 
with respect to the additive quark mass renormalization 



A/,, 



(w)/,o-f (w)s,o' 



(5.1) 



We normalize this combination by its zero temperature value to cancel the multiplicative renormal- 
ization factor. In Fig. H we show on the left hand side results for niq = O.lm^ with A^^ = 6,8 as 
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Figure 4: 

tibility is 
and nig = 



We plot Ai s for the two different light quark masses (left). The disconnected light chiral suscep- 
shown on the right hand side for light quark masses niq = Q.lnis (circles), niq = O.lmj (squares), 
0.05m.s (triangles). 



well as results for = 0.05mv with A^^- 
for both data sets with m„ 



= 8. We observe a sharp drop of A/^^ in the region of T^. 
O.lnis and = O.OSm,. We note that we will present data for more 
different light quark masses in the future which will allow us to probe the limit 0. With only 
two different light quark masses we cannot investigate this limit at the moment since near Tc terms 
proportional to rriq as well as ^/rn^ are present. 

On the right hand side of Fig. ^ we plotted the disconnected light chiral susceptibility Xi.disc 
for different light quark masses. We find that Xi.disc shows a strong quark mass dependence over 
a wide temperature range. This can be understood as an increase in the fluctuations of Goldstone 
modes below T^. The transition temperature should then be found near the right edge of the peak. 
Further evidence for this picture is given by the observation that the height of the peak of Xi,disc 
scales as 1 / ^JrPi^. 
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Summary 

We presented new data for the equation of state for niq = 0.05ms close to the transition tem- 
perature. The pion mass is as low as « 160 MeV. We compared the trace anomaly to earlier results 
for ruq = O.lmj. A total shift of the transition region of about 4 MeV is found towards smaller 
temperatures. About 0.5 MeV of that shift can be contributed to the shift in the scale. Compar- 
ing to the Hadron Resonance Gas model (HRG) we find that at larger temperatures the results for 
niq = 0.05ms are consistent with HRG model while there are deviations at lower temperatures. 
Calculations at different temperatures are currently ongoing. 

We have studied the behavior of A/^^ for the two different quark masses. A sharp drop in the 
vicinity of the critical temperature is observed. Furthermore, we have presented results for the 
disconnected Ught chiral susceptibiUty. In the transition region Xi,disc shows a broad peak for the 
different light quark masses which rises as 1 1 ^Jm'q. This gives support for a picture where the 
Goldstone modes are causing the fluctuations below Tc- 
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